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^ Abstract 

We develop a theory of adiabatic response for open systems gov- 
erned by Lindblad evolutions. The theory determines the dependence 
of the response coefficients on the dephasing rates and allows for resid- 

1-1 ual dissipation even when the ground state is protected by a spectral 

gap. We give quantum response a geometric interpretation in terms of 
£> Hilbert space projections: For a two level system and, more generally, 

for systems with suitable functional form of the dephasing, the dissi- 
pative and non-dissipative parts of the response are linked to a metric 
and to a symplectic form. The metric is the Fubini-Study metric and 

00 the symplectic form is the adiabatic curvature. When the metric and 

symplectic structures are compatible the non-dissipative part of the 
inverse matrix of response coefficients turns out to be immune to de- 
phasing. We give three examples of physical systems whose quantum 

• »-j states induce compatible metric and symplectic structures on control 

rN space: The qubit, coherent states and a model of the integer quantum 

& Hall effect. 

Two frameworks that provide insight and understanding of transport coef- 
ficients are Kubo's theory of linear response [HE] and the theory of adiabatic 
response [3], IH E]. Both have a quantum version and a classical version, agree 
when there is overlap, and endow the non-dissipative transport coefficients 
with the geometric meaning of the adiabatic (Berry's) curvature [6]. 

A notable success of this approach to "geometrization of quantum re- 
sponse" has been its application to the integer quantum Hall effect where 
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one observes quantized resistances of the form h/ne 2 with n an integer. The 
accuracy and robustness of the quantization is understood as reflecting the 
nontrivial topology of the quantum state characterized by a topological in- 
variant, the Chern number, which is the n measured as the Hall resistance 
[7] . "Geometrization of transport" also lies at the heart of topological clas- 
sification of states of matter [HI |TU] and the considerable current interest 
in topological insulators [IT] . 

Our aim here is to carry over the program of "geometrization of quantum 
response" to open quantum systems. More precisely, extend the theory of 
adiabatic quantum response from the Hamiltonian setting to the Lindbla- 
dian setting [T2J H3]. The extension gives geometric meaning to both the 
non-dissipative and dissipative response coefficients, and allows to examine 
how they are affected by decoherence. It allows to address the stability of 
topological quantum numbers, such as Chern numbers, against dephasing. 

In adiabatic response, the Hamiltonian, H(<fr), is viewed as a function of 
control parameters <fi which drive the system [31 HI [T5] . The space of controls 
shall be denoted by hA. We focus on observables which are gradients of H 
over Ai, namely 

For example, in the system shown in Fig. [TJ the controls are magnetic fluxes 
which by the Aharonov-Bohm effect may be thought of as angular variables 
so that control space M. is the torus T 2 . is the loop-current in the /i loop 
as one can see by the principle of virtual work. For other notions of currents 
in the Lindbladian context, see [To] . 

Our main object of interest is the matrix of response coefficients / defined 
through the linear and instantaneous terms in (f> in the (adiabatic) expansion 
of the response 

Ti(p t F,) = f, v ^{t) + ... . (2) 

Summation convention over repeated indices is implied. All other terms, 
which are not proportional to 0, do not concern us. In the example of Fig. [l] 
cf)^ is the emf on the v loop and / is then the conductance matrix relating 
loop currents to emf's. 

The quantum state pt in Eq. ^ arises from the ground state as the 
solution at time t of the equation of motion 
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Figure 1: A torus in coordinate space threaded by two loops (red) carrying 
fluxes <pi and 02 viewed as a model for the QHE. The (strong) magnetic field 
(not shown) is perpendicular to the surface of the torus. The control space 
Ai is the space of fluxes. By the Aharonov-Bohm periodicity, Ai is a torus 
as well. 

governed by the adiabatically changing Lindbladian: 

c( P ) = -i[H, P ] + ^T([r>,r;] + [r a;P r*]) ; (3) 

a 

T a = corresponds to a unitary evolution. 

In adiabatic transport [3J HI [3] one is interested in the situation where the 
control parameters <j) move adiabatically along a path in control space Ai. It 
turns out that adiabatic methods used to study the unitary case, r a = 0, can 
be extended to also study open systems described by Lindbladians, provided 
the instantaneous stationary states move continuously with the controls. The 
response of a stationary state, C(p(4>)) = 0, of the Lindbladian £ to a driving 
is given by (see Appendix B) 

f, v = Ti{F,C-\d uP )). (4) 

/ M „ has a geometric interpretation when the stationary state p is a (spec- 
tral) projection. A particular choice of Lindbladians that achieve this is 
r Q = T a (H) for some function of H. We call this family of changing Lind- 
bladians, where T a (H) are slaved to the instantaneous Hamiltonian, dephas- 
ing Lindbladians. The family describes processes that conserve energy and 
entropy increases without heat exchange. In appendix A we give an example 
showing how dephasing Lindbladians naturally emerge in certain stochastic 
unitary evolutions. 
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The slaving of the decoherence terms to the instantaneous Hamiltonian 
has the consequence that the instantaneous stationary states of the Hamil- 
tonian and the Lindbladian coincide. Moreover, when the controls (f) vary 
adiabatically, one expects spectral projections of H to evolve so that they 
remain close to the corresponding instantaneous spectral projections P{4>)- 
The evolution of spectral projections can be described in geometric terms. 
Geometric quantum response is concerned with the relation between the re- 
sponse matrix / and the geometry of P{4>) on control space. 

Before we describe this relation, it is convenient to review the geometry on 
the control space associated with a single spectral bundle \-> P((f>). Con- 
sider the (operator-valued) 1-form giving the natural adiabatic connection 
(here P± = 1 - P) 

A = A fl d<f>>*, A ll = P ± d ft P (5) 
from which one can construct the second rank tensor on control space 

2 Ti(A ® A*) = ( gtlv (4>) - i uv(0))tty" ® (6) 

where <g> is the product of forms. The symmetric part is the natural notion 
of infinitesimal distance for projections, the Fubini-Study metric [TT] . 

9,u = Tr P ± {d u P, d,P} = Tr (d u P) (d,P). (7) 

The antisymmetric part gives the adiabatic (or Berry's) curvature [6], 

uv = iTr(P ± [d tl P,d u P}) = -tTr(P[d,P,d u P]). (8) 

The basic properties of g and cu are direct consequences of the Eqs. (]7| [8p 
and are summarized in the following statement. 

Proposition 1. g > defines a metric on the space of controls Ai. It is 
even under time-reversal and it is also even under electron-hole exchange, 
g{P) = g(P±). u is a closed 2- form, doj = 0, endowing M. with a symplectic 
structure, if non-degenerate. It is anti-symmetric under time reversal and 
hence u = if H{4>) is time reversal invariant, u is anti- symmetric under 
electron-hole exchange, oo(P) = —oj{P±). 
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Geometry of quantum response: 

For the sake of simplicity we consider the case where H is finite dimen- 
sional matrix and T a (H) are real functions of H. The (instantaneous) eigen- 
values of H shall be denoted by Ej, and the corresponding eigenstates by 
The ground state shall be denoted by j = 0. The eigenvalues of C 
(corresponding to eigenstates \j)(k\) shall be denoted by A^. By inspection 
Im(Ajfc) = — £j + Ek while Re(Ajfc) < 0, i.e. all eigenvalues lie in the (closed) 
left half-plane. If two states are degenerate Ej = Ek then clearly \jk = 0. 

It will be convenient to introduce the notion of dimensionless dephasing 
rates 7^ > associated with the pair of non-degenerate eigenstates 

7* = -^4 > 0. (9) 
\£j ~ £k\ 

For a pair of degenerate states, where both the numerator and denominator 
vanish, the dephasing rate is defined to be zero. 

We are now at a position where we can state our first main result. In 
Appendix B we describe a formula relating the response matrix / to the 
spectral projections and their gradients for general dephasing Lindbladians. 
The formula admits a particularly simple geometric interpretation provided 
we make a specific choice for T(H). We shall first state the result in this 
special case and then comment on the general case. 



Theorem 1. The matrix of transport coefficients f , Eq. associated to 
the adiabatic evolution characterized by a single dephasing rate, Eq. |s|) ; 7 = 
Ijo > to all other energy levels j , is given by 

1 + 7^ 1 + / y z 

g is the Fubini-Study metric of the ground state bundle and u> its adiabatic 
curvature. 

Let us make the following observations: 

1. The condition of a single dephasing rate is automatically satisfied in any 
two level system. It is also satisfied if H has two degenerate eigenvalues. 
When H has 3 or more distinct eigenvalues the condition is non-trivial 
and can be interpreted as a condition on the functional form of T(H), 
e.g. T(H) = y/j(H — Eq) with e the lowest eigenvalue of H. 
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2. The dissipative response is associated with the symmetric part of / and 
is fully determined by the metric g. It vanishes in the limit 7 = as 
it must — there is no dissipation in adiabatic unitary evolution with a 
gap condition. In contrast, in open dephasing systems, a gap condition 
does not provide protection from dissipation. 

3. The non-dissipative response is associated with the anti-symmetric part 
of / and is fully determined by the adiabatic curvature and the dephas- 
ing rate. (An extension of the Berry phase to systems with decoherence 
in the form of a complex phase is studied in [13] . ) 

4. The theorem gives a geometric interpretation to both the dissipative 
and non-dissipative parts of the adiabatic quantum response. 

5. For weakly dephasing systems, < 7 >C 1, the dissipative response de- 
pends linearly on 7 while the non-dissipative response depends quadrat- 
ically on 7. Dephasing affects both the dissipative and non-dissipative 
response coefficients. 

6. Chern numbers are integers obtained by integrating the adiabatic cur- 
vature over a closed 2-dimensional control space. The theorem says 
that the relation between the control space average of transport coef- 
ficients and Chern numbers is a function of the dephasing 7. Hence, 
topological quantum numbers are not robust against dephasing. 

7. The theorem generalizes to the multi-dephasing rate case, where 7-,o of 
Eq. (|)J is j-dependent, at the price of replacing P±d fl P by a weighted 
sum of Pjd^P. The Fubini-Study metric is then replaced by a metric 
that does not have a standard name. 

8. The linear response formula, Eq. ([2]), gives the leading term in the 
adiabatic expansion of the response provided the spectrum of H(<f>) 
is independent of (p. In the general case, when the eigenvalues are <ft 
dependent, the expansion Eq. ^ has additional terms at low orders: 
One which is not small, but depends on only, and one which is of 
the same order as (pit), but is given by a quadratic expression in cf)(t') 
integrated over t' < t. However, as these terms are not proportional 
to the instantaneous driving (fi(t) we shall not consider them here (cf. 

i)- 
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Control space with compatible metric and symplectic structures: 

Our second main result concerns the special class of Hamiltonians when 
the metric g and the symplectic structure u are compatible. In the case of 
a single pair of controls this is expressed by det g = det uj. In general and 
in terms of a basis in which both g and u are 2 x 2-block diagonal, the two 
structures are compatible if they are so inside each block. Equivalently, we 
say that g and uj are compatible if 

uj- 1 g + g- 1 oo = 0. (11) 



Eq. (11) implies the equality of the determinants. 

By Proposition [TJ compatibility is possible only when time reversal is 
broken. Below, we shall give three natural physical examples which give 
rise to compatible metric and symplectic structure. The three examples 
correspond to the three prototypical control spaces M.: the sphere, the plane 
and the torus, respectively. 

Theorem 2. Suppose that g and uj are compatible. Then the inverse of the 
matrix of quantum response of Theorem [7] has the form 

f-^-yg^+u- 1 . (12) 



The claim is easily verified by multiplying Eq. (10) by Eq. (12) and using 



Eq. (11). 



Remarkably, the non-dissipative response associated with the antisym- 
metric part of f' 1 is independent of 7 and so immune to dephasing. It is 
determined by the adiabatic curvature alone, just like in the case of unitary 
evolution. 

Compatibility test: 

Given P one is interested in simple tests that tell whether g and uj are 
compatible without explicitly computing the matrices g and uj. 

Theorem 3. The following are equivalent. 



(i) g, uj defined on j\A. are compatible in the sense of Eq. (11). 



(ii) M. has local holomorphic and antiholomorphic coordinates z^ , z^ mak- 
ing it into a complex manifold and the map P(4>) satisfies P±djP = 0. 
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(Hi) The image of the map P : Ai — >■ Qi(n + l,r;C) a complex sub- 
manifold of the complex Grassmannian manifold Gr(n + 1, r; C) . (P/ere 
r = rank(P) and n + 1 is the dimension of the Hilbert space.) 



Moreover when P is a one dimensional projection (i.e. r = 1) then (ii) is 

WW 



equivalent to the claim that P may be expressed as P = Tyny where 



is holomorphic, i.e. dj \ip) — 0. 

We show how the test (ii) implies compatibility (i) when Ai is two di- 
mensional. We write locally z = cf) 1 + r<f) 2 , where r = T\ + ir 2 with r 2 7^ 0. 
Then B = 2^(t<9i — d 2 ) and the assertion (ii) is equivalent to the assertion 
that 

tA x - A 2 = 0. (13) 
Substituting this identity in Eq. ^ gives 

922 = \r\ 2 gu, 9x2 = ngu, u 12 = T 2 g u - (14) 

The equality of the determinants now follows by inspection. 

For P = |0) (0|, (ii) is equivalent to d s \0) oc |0) and hence the possible z 
dependence of |0) 2 is merely through a normalization/phase factor. 

Complex structure-Kahler manifolds: 

The compatibility condition is equivalent to the statement that the op- 
erator J = u~ 1 g acting on T^Ai, the tangent space at <ft, satisfies J 2 = —1. 
It can therefore be used to turn the real vector space T^Ai into a complex 
vector space by defining multiplication by the scalar a + i(3 G C to be given 
by the action of the the operator a + J. 



The compatibility condition Eq. (11) automatically holds for the Fubini- 



Study metric and the adiabatic curvature when the control space is the whole 
projective space PC n and the projections P parametrize themselves. This is 
related to the fact that PC™ is a Kahler manifold. The same also holds for 
the Grassmanian manifold. Our control space Ai may be viewed (through 
the map <p ^ -P(0)) as immersed in PC™ when rank(P) = 1 or in the Grass- 
manian when rank(P) > 1. Indeed we defined g M and u M on Ai using this 
identification by restricting ^ PC " ;C( ; IPCn to its immersion. g M ,u M inherit the 
compatibility from g " , uj vcn if and only if Ai is a complex submanifold of 
PC™ (or of Gr(n + 1, r)), in which case Ai inherits the Kahler property, too. 
This leads to the Theorem [H 
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We conclude by giving three examples of natural controlled physical sys- 
tems whose ground state bundle has a metric compatible with the curvature. 

The Qubit family: The controlled Hamiltonian of a qubit is 

H = (f>-a, 0G§ 2 . (15) 

The control space Ai = § 2 is the unit sphere. The ground state projections 
are given by 

2P = l-0-(T. (16) 

The associated Fubini-Study metric and the symplectic form are readily com- 
puted from Eqs. (pTT [8l 



2<7ij = di4> ■ dj(f), 2uij = —0 • di<j) x dj4>. (17) 

Both give consistent areas (half the standard area) on § 2 = At. 

Coherent states: The Hamiltonian of a phase-space controlled oscillator 

is 

H{C^) = \{v-tf + \{x-Q 2 (18) 

where ((, fi) G 1R 2 = Ai, the Euclidean plane. The manifold of ground states 
is the coherent states. The controls are boosts and shifts hence 

A x = P±d ( P = -iP±\p, P] = -iP±{p - n)P. (19) 

Similarly, 

A 2 = P ± d tM P = iP ± [x,P] = iP ± (x-()P (20) 

Since a = (x — () + i(p — //) annihilates the ground state of the shifted 
oscillator, Eq. (13) holds with r = i 

P ±a P = 0. (21) 

Landau Hamiltonian: As a third example of a control space with Kahler 
structure consider the Landau Hamiltonian associated with a torus of unit 
area and skewness r = n + ir 2 , r 2 > 0; (see Fig. [2]). The torus is threaded 
by Aharonov-Bohm magnetic fluxes (0i, <\) 2 ) (see Fig. fl| and is penetrated by 
constant magnetic field 2nB, B an integer, (see Fig [2J) . The corresponding 
Hamiltonian is [18] 

H(t,4>) = \d*D (22) 

T 2 
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Figure 2: A flat torus with skewness r is obtained by identifying the parallel 
ends of the skewed parallelogram. A constant magnetic field with an integer 
number, B, of flux quanta penetrates the torus. 



where D = i(rd x — d y ) — 2nr(By + 0) and cf> = (f>i — <fe/T. We impose the 
usual magnetic translation boundary conditions [19|: 

i/;(x, y) = ^{x + l,y) = e 2mBx <i[>(x, y + 1). (23) 

The (single particle) ground state is 5-fold degenerate with energy E = 
independent of r and 0. For simplicity we consider the case of a single particle 
and B — 1. The lowest Landau level is then related to theta functions [20J: 

oo 

ip(x,y) = (2/r 2 ) 1 /V" |r|2(Im * )2/r2 e 2winx e inT{v+n+ ^ )2 . (24) 

n=— oo 

Since, except for the normalization, ip depends only on the "holomorphic 
coordinate" and not on the "anti-holomorphic" coordinate (f) it follows 
from Theorem [3] that g and u are compatible. Alternatively, the analyticity 
of ip is also seen from that of D; note that the boundary conditions are 
independent of <fi. 

Conclusions: 

This work studies the response of driven open quantum system governed 
by adiabatically evolving Lindbladians. When the Lindblad operator is de- 
phasing, coherence and phase information is degraded but the energy of the 
system is still conserved. In this case, we find that the response admits a geo- 
metric interpretation induced from the behavior of the instantaneous Hilbert 
space projections. 

We focus on the response associated with observables that can be derived 
by the principle of virtual work from the Hamiltonian of the system. For these 
observables we find that, in contrast with the case of unitary evolution, there 
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is residual dissipation: A spectral gap does not protect against dissipation in 
open systems in general. 

Our first main result, Theorem [TJ concerns the special case of dephasing 
two-level systems and, more generally, certain dephasing systems character- 
ized by a single rate 7. In this case, the response matrix admits a simple 
geometric interpretation: The non-dissipative transport is proportional to 
the adiabatic curvature uj and the dissipative response to the Fubini-Study 
metric g. Since the proportionality factor relating non-dissipative transport 
to the adiabatic curvature is a function of the dephasing rate 7, it follows 
that the robustness of Chern numbers does not translate to a robustness of 
the transport coefficients against dephasing. 

Our second main result, Theorem |2j is concerned with a geometric mecha- 
nism which provides protection against dephasing for certain non-dissipative 
response coefficients. The mechanism leading to such a protection is the 
compatibility of the adiabatic curvature with the metric. In this case, the 
relation between transport coefficients and Chern numbers is independent of 
the dephasing rate. 

Our third main result concerns test of compatibility. In particular, we 
show that when the instantaneous stationary state are holomorphic functions 
of the driving parameters, the Fubini-Study metric is compatible with the 
adiabatic curvature. 

We conclude with three fundamental examples of compatible physical 
Hamiltonians associated with the three types of control spaces. The first 
example of a compatible system is the qubit Hamiltonian whose control space 
is the sphere, the second example of a compatible system is the Harmonic 
oscillator whose control space is the plane and the third example is a model 
of the quantum Hall effect on a torus. 

Acknowledgments. We thank Misha Reznikov for suggesting that the 
Hall resistance is better behaved than the Hall conductance. This work is 
supported by the ISF and the fund for Promotion of research at the Technion. 
M. F. was partially supported by UNESCO fund. 

Appendix A 

Here we describe an example illustrating how adiabatic dephasing Lindbla- 
dian, with slaved dephasing term, naturally arise from stochastic unitary 
evolutions. The example is a stochastic variant of Berry's paradigm of the 
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notion of adiabatic curvature, namely, a spin 1/2 in a magnetic field [5]. The 
evolution equation is 

p = -i[B-a,p] (25) 

where B G IR 3 is a time dependent magnetic field and a the vector of Pauli 
matrices. The case considered by Berry is when B changes its orientation adi- 
abatically say with fixed magnitude. We want now to consider the stochastic 
version of this model where the magnitude of B is a stochastic variable while 
its orientation is changing smoothly (adiabatically) in time. Formally, this 
corresponds to replacing B in the evolution equation by 

B -> W t B 

where W t is (scalar, biased) white noise. The canonical interpretation of 



Eq. (25) as a stochastic differential equation goes through the Ito calculus 
|21j . To do so, it is convenient to expresses white noise in terms of the 
corresponding Brownian motion 

db t := W t dt. 

The rules of Ito calculus say that dp has to be expanded to first order in dt 
and to second order in db. This gives the stochastic evolution equation 

dp = -t[H ,p]db - l(db) 2 [H , [H ,p]], H = a-B 

where B is the smooth (non stochastic) function of time. In particular, 
it follows that the (noise average) state p a = E(p) satisfies the adiabatic 
Lindblad equation 

p a = C{p a ) = -ip[H ,p a ] - \D[Hq, [H 0l p a ]] 

where p is the bias of the white noise p = E(W t ) and D its variance 
E(WiWs) = DS(t — s). If D ^ 0, this gives a dephasing evolution where 
the dephasing is slaved to the time dependence of the Hamiltonian. (A gen- 
eral framework for deriving Lindbladian for general stochastic evolutions is 
described e.g. in 



Appendix B 

Here we outline the proof of Theorem [T] by evaluating the terms proportional 
to in Eq. (B. 
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Let Pj denote the spectral projections for H. Since C(Pj) = 0, the 
spectral projections are instantaneous stationary states. Let P = Po denote 
the projection on the ground state. We also denote Ejk = \j) (k\. This is an 
eigenvector of the Lindbladian with eigenvalue Xjk i.e. C(Ejk) = XjkEj^. 

By the adiabatic theorem the states adheres to the spectral projection, 
p(t) = P((p(t)) + 0(<f)). The first order correction Sp to the state satisfies 

C(5p) = P, (26) 

as can be seen from the substitution p = P + 5p into the Lindblad Eq. ^ and 
using C(P) = 0. The correction can be decomposed as 5p = 5±p + 5\\p into 
parts 5±p G Range £ and <5||p e Ker£, which are orthogonal with respect 
to the inner product defined by the trace. Note that £ considered as a 



map on Range £ is invertible and that P G Range £. Thus Eq. (26) implies 
5±p = £ _1 (P), where the inverse C~ l is well defined. In fact, since the 
eigenstates of £ are Ejk, one may readily write 

5 ±P = C-\P) = J2 {jl P^E jk . (27) 

Strictly speaking we restricted here to the case of simple eigenvalues; more 
generally (j| P \k) — between degenerate eigenstates, whence the appropri- 



ate reading of the sum (27) is by omitting such pairs. 

The complementary part 5\\p may be determined as well (cf. jl]) and 
happens to depend on history, but will not be needed. 

Now, p carries two contributions to the response Eq. ([2]), of which the 
leading one, Tr(PF M ) equals d^Ti{PH) = d^eo. This term is not propotional 
to and does not concern us (note that it vanishes when the spectrum is 
independent of 0, cf. observation 8 after the theorem). As for the first 
order correction 5p, two contributions arise in turn through F, 



dH 



f 



d^^SjPj. The first, YK®v £ j)Pji ^ es m Ker£ and matches 5\\p. This term 
does not concern us either and again vanishes when d^Ej does. The other 
part ^2 £jdfj,Pj lies in Range £ and gives the requisite linear response term of 
the expectation value (F^) = Ti(pd^H): 

^^((fyPOM =J2^i-^)Tr((d fl P i )S ± p), (28) 

i 
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where we used <9 M J2i Pi = 0- Eq. (28) can now be written as J2i( £ i ~ £ o)Ai, 
where 

A = £(fc|^|j}^. (29) 

Using 

01 P \k) = (j\ PjPP k \k) = (4,o + S jl0 ) (j\ P \k) (30) 



it follows that the double sum in Eq. (29) reduces to the single sum 

A=E (MMMM +K ), (si) 

Since Xkj = \jk we get A% is manifestly real as it must be. Using the fact 
that (recall that i ^ 0) 

(0| d^Pi \j) = (0| Pd,P \j) = - (0| (d,P)P t \j) = -Sij (0| d,P \j) (32) 

we finally find 

J> - e ) A' = E Tr {P{d li p)P j P) + c.c. 

where 7^0 > is the dimensionless characterization of the spectral data of 
Eq. (§. 

Simplification occurs for jjo is independent of j. This is, of course, au- 
tomatically the case for a two level system where j takes one value j = 1. 
(Similar simplification occurs when there is one dominant The sum 

over j can now be carried out explicitly 

(F,) = - £oH- = ^{(d,P)P±P) + c.c. (34) 

Writing P = ^2 u (d u P) <f> we obtain the expression in the theorem. 
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